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Abstract. We show how the buildup of halos by merging satellites 
forces an inner cusp, with a density profile p oc r~" where a — > Oa ~ 1- 
Our analysis is based on a new prescription for tidal stripping as a func- 
tion of a{r), using a simple toy model which matches N-body simulations. 
In a core of a < 1 there is tidal compression rather than stripping and 
the satellites sink towards the halo center, causing a rapid steepening of 
the profile to a > 1. Where a > 1, the stripping of each satellite shell is 
preceded by gradual puffing up, which makes the stripping more efficient 
at larger a, causing flattening where a is large enough. Therefore, we 
can show using linear perturbation analysis that a sequence of mergers 
slowly leads to a fixed point a{r) = Oa- This result implies that a cusp 
is enforced as long as enough satellite material makes it into the inner 
halo and is deposited there. We conclude that in order to maintain a flat 
core, as indicated by observations, satellites must be disrupted outside 
the core, e.g., because of puffing up due to baryonic feedback effects. 



1. Introduction 

The 'standard' model of cosmology, CDM, which assumes hierarchical buildup 
of structure, is facing difficulties in explaining observed properties of galaxies, 
such as the number density of dwarfs (e.g., Klypin et al. 1999b), the angular- 
momentum crisis (e.g., Navarro & Steinmetz 2000), and the cusp/core prob- 
lem. Our approach in addressing these problems within CDM is to isolate and 
model in simple physical terms the key relevant processes, as a guide for possi- 
ble solutions. We first model the buildup of dark-matter (DM) halos in N-body 
simulations based on tidal effects, and then incorporate the inevitable baryonic 
feedback processes in an attempt to explain the apparent discrepancies. We 
address the angular momentum problem in Mailer & Dekel (2002) and Dekel 
& Mailer (2002), and summarize here our progress in the cusp/core problem 
(Dekel & Devor 2002; Dekel et al. 2002). 

Cosmological N-body simulations have revealed that the density profiles of 
DM halos scatter about a universal shape, p{r) = (?^/rc)~" (1 + r/rc)"~^, 
with an inner cusp of slope —a. Navarro, Prenk & White (1995, NFW) found 
this function, with a ~ 1, to be a good fit in the range (0.01 — 1)-Rvir for 
different hierarchical cosmological scenarios. High-resolution simulations of a 
few individual halos (Moore et al. 2001; Klypin et al. 2001) found that the cusp 
could be as steep as a ~ 1.5, though it may flatten towards a ~ 1 at r < 0.01i?vir 
(private comm. with Navarro, Frenk, Springel & White). While the formation 
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of a cusp with 1 < a < 1.5 has been estabUshed in the simulations, a basic 
understanding of its origin is stiU lacking. An even more intriguing puzzle is 
introduced by observations of low surface-brightness galaxies, whose centers are 
dominated by their DM halos, which indicate flatter inner cores with a ~ (de 
Block et al. 2001). This seems to introduce a challenge to the CDM paradigm. 

In the first two sections we develop a toy model for tidal stripping and test 
it against an N-body simulation. In §2 we describe the compression and rapid 
steepening in a core. In §3 wc address stripping where a > 1. In §4 we analyze 
the convergence to an asymptotic profile. In §5 we discuss our results. 



2. Core into Cusp due to Tidal Compression 

We consider a fixed spherical halo of mean density profile p{r) oc M(r)/r^. De- 
note a{r) = —din p/dlnr such that locally p oc r~", with a{r) constant or 
increasing in the range < a < 3. Consider a satellite at r, spiraling into the 
halo under gravity and dynamical friction. The maximum tidal force by the halo 
on a unit mass at satellite radius £ <C r) is 

Ftide = ^^, ^i{r) = 2M{r)-r— = {a-l)M{r). (1) 

The maximum is obtained along the line connecting the halo centers, but while 
the radial component is smaller in all other directions, it is of the same sign. 
The familiar case is of a point mass, a = 3, for which the tidal force is pulling 
outwards from the satellite. For flatter halo slopes it becomes weaker oc (a — 1). 
An important feature that is often overlooked is that when the halo density 
profile is flat enough, a < 1, the tidal force reverses direction into compression, 
resulting in accretion rather than stripping. Note that the critical slope of unity 
coincides with the cusp of the NFW profile, and that in a core, a ~ 0, the tides 
induce strong compression. 

This effect can be demonstrated using an N-body simulation of a merger 
(following Mihos k, Hernquist 1996). We use a total of 10^ DM particles, the 
mass ratio is 1:10, the satellite spirals in on a quasi-circular orbit, and the initial 
profile is a truncated isothermal sphere with a core, where a ranges in practice 
from 0.6 to 2.8. In Fig. 1 (left), the stripping of each satellite shell, marked 
by the onset of a steep rise, can be identified with a halo radius r, with a{r) 
decreasing from 3 downwards as one moves from outer shells inwards. However, 
as the stripping point approaches a ~ 1, inwards to the 30% satellite mass shell, 
the stripping stops and some shrinking can be seen instead. 

The compression at a < 1 implies that any part of the satellite which 
makes it intact into the halo core would sink towards the center without further 
stripping. This should cause a rapid steepening of the profile to a > 1, as seen 
in Fig. 1 (right). The NFW inner slope of a = 1 is thus a robust lower bound, 
as we know from cosmological simulations; a flatter density core cannot survive 
as long as satellites deposit enough mass in the inner halo. 
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Figure 1. Left: Evolution of the satellite mass profile during the 
merger simulation. Shown are the radii of spheres about the satellite 
center, encompassing given fractions of the satellite mass. Marked 
along the x-axis are values of a(r) corresponding to the momentary 
position of the satellite center in the halo. The inner 25% of the satellite 
mass is never torn away — it enters intact the a < 1 core. Right: 
Halo density profile before (solid) and after (dashed) the merger. The 
compact inner satellite sinks into the halo core and produces a cusp 
with Q > 1. 



3. Tidal Pufiing-Up and Stripping at a > 1 

At a > 1, the effects are more subtle. Let £ mark shell radii within the un- 
perturbed satellite, whose mean-density profile is a{£) oc m{£)/t^ . Assume that 
when it is at halo radius r, mass is lost beyond a momentary stripping radius 
i{r) and is added to the halo at r (on average). We wish to determine the 
correspondence between I and r at stripping. Define ip{r,i) = p{r)/a{i). Tra- 
ditionally, the stripping radius is assumed to be determined by the resonance 
condition ■ilj{r,i) = 1, but this ignores the earlier effects of tides on the satellite 
structure. A key new feature in our analysis is that, as r decreases, the tides 
stretch the satellite orbits, which can be modeled as an effective puffing up of the 
relevant shells before they are being torn away. We define for every shell i when 
the satellite is at r a momentary puffing factor by p{r,i) = £p/i, where ip{r) is 
the momentary shell radius. One can then show that the resonance condition, 
for r and i at stripping, becomes 

= a'^p-^ . (2) 

In the regime where a > 1, we expect puffing, p > 1, so the corrected resonance 
condition implies ip < 1, and for large a even -0 <C 1. This means more efficient 
stripping compared to the old condition ignoring puffing. To obtain an explicit 
stripping condition we wish to express ip asa function of a, so we need to estimate 
how p{i) evolves as the satellite falls into smaller r positions. By applying an 
adiabatic invariant, we obtain a puffing equation for any shell i when it is at r: 

p- {a-l)ipp'^ = 1. (3) 

In the two equations above we have omitted geometrical factors of order unity, 
to be calibrated later using simulations (see Dekel & Devor 2002). When we 
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Figure 2. Left: Density ratio V at stripping as a function of halo 
slope a. The predictions of the model are shown (dashed) for three 
different choices of the geometrical factors. The no-puffing prediction 
is '0 = 1- The result of the N-body simulation (solid) demonstrates 
that the puffing model can provide a good approximation. Right: A 
schematic illustration of satellite mass deposit in the halo. Shown are 
an NFW halo profile p{r), and a homologous satellite a{i) properly 
shifted to the left and upwards. The arrows connect shell radii i to the 
halo radii where they are deposited r. The horizontal dashed arrows re- 
fer to stripping when puffing is ignored, ^{a) = 1. This would steepen 
the proffie, as steep regions of a{i) are deposited at flatter regions of 
p{r). The solid arrows illustrate realistic stripping after puffing. The 
vertical displacements, which grow with r, refer to ip{a) < 1. The slope 
at i may be flatter than the slope at r such that the mass tends to be 
deposited at larger r and the result is flattening of the halo profile. 



combine the above equations we obtain a new stripping condition: 

na[r)\ jo.l/(a-l) a > a, ' 

Given the slope profile a(r), it relates every satellite shell £ to the position r 
where it should be stripped. The function ip{a) is thus predicted to decrease 
monotonically as a function of a, towards values of order 0.1-0.2 at a ~ 2 and 
below 0.1 at a = 3, Fig. 2 (left). This implies that the stripping process is more 
efficient for steeper halo profiles. We also implement in the deposit prescription 
additional effects near and below a = 1, due to the finite size of the satellite 
and the tidal compression in the inner halo. These effects may cause deposit of 
satellite material before the stripping condition is fulfilled. 

The puffing before stripping can be seen in the merger simulation, e.g., for 
the outer and intermediate shells in Fig. 1 (left). The magnitude of the puffing is 
30-50%, as expected, corresponding to a factor of 2-3 in density. A similar effect 
has been qualitatively noticed in simulations before (e.g. Klypin et al. 1999a, 
Fig. 6). In the merger simulation, we measure the deposit radius r(i) from 
the final distribution of stripped satellite mass about the halo center [can be 
read from Fig. ?? Fig. 1 (left)]. The corresponding values of ip{r,tj are plotted 
against a{r) in Fig. 2 (left). The qualitative agreement between the simulation 
result and the model predictions indicates that despite the crude approximations 
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made, our very simplified model mimics the main features of the tidal stripping 
and deposit process. When puffing is ignored, V' = !> the model clearly fails. 

4. Halo Asymptotic Profile 

If the stripping is described by a condition similar to eq. ??, with ip{a) a de- 
creasing function, the profile evolves slowly towards an asymptotic stable power 
law with aa ~ 1- We assume that the halo and satellite are drawn from a cos- 
mological distribution; they are homologous, with their characteristic radii and 
densities scaling like ic/rc oc m^^^'^^^^ and Cc/Pc oc m~'^ , where ly ~ 0.33 for 
ACDM. Fig. 2 (right) helps understanding the origin of an asymptotic slope. 
We write p{ina\{r) = p{r) + a{£)£^ /r^, and obtain for the change of a in a merger 

One can see that every power law is a self-similar solution, Aa(r) = 0, but not 
necessarily a stable one. When a is increasing with r, i/r is decreasing with 
r. Thus, when puffing is ignored, tjj =const., one has continuous steepening, 
Aa{r) > 0. With realistic puffing, l/ip is increasing with r, which can produce 
a stable fixed point at a certain asymptotic value where Aa = and the 
second derivative is negative. A rigorous linear perturbation analysis determines 
the rate of convergence to Oa and yields an equation for its value for a sequence 
of mergers with the same mass ratio: 

Aa cx a{a - 3)ip' {a)/ip{a) + 3ln[{m/Myi;{a)] = . 

The solutions are typically in the range 1 < Oa < 1-5. 

In order to test the linear analysis, we perform toy simulations of the profile 
buildup by mergers, where we implement the stripping recipe of §3 (replacing 
the crude stripping recipe used in earlier work, e.g., Syrc & White 1998). For 
given halo and satellite profiles, we solve the stripping equation for i{r) and 
add the stripped satellite mass to the halo accordingly. Among other numerical 
complications, wc implement a smoothing scheme to ensure that a(r) remains 
monotonic. Wc then follow a sequence of cosmological mergers and study the 
evolution towards an asymptotic slope. Fig. 3 (left) shows the convergence of a 
at a fixed r to the asymptotic value. Fig. 3 (right) shows how the profile evolves 
through momentary profiles which are probably more relevant for comparison 
with real halos at different times during their buildup process. These figures are 
for certain given values of the geometrical factors, the mass ratio and u. In Dekel 
et al. (2002) we address a sequence of mergers with a cosmological distribution 
of mass ratios, and the robustness to the cosmological model. 

5. Discussion 

Our analysis demonstrates that the way to maintain a flat core is by disrupting 

satellites outside the core. This may be achieved if the cores of satellite halos 
are puffed up due to gas processes. As an example, the following speculative 
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Figure 3. Left: Toy-simulation evolution of slope a at r = O.lrc due 
to a sequence of mergers n = 1,600 with mass ratio m/M = 0.3. The 
initial profile is a generalized NFW with a either zero or 2. When a < 
1, the slope steepens rapidly to a > 1 within a few mergers (§2), and 
then it converges slowly from either side towards an asymptotic value 
(§4). Right: Corresponding evolution of slope profile a{r), starting 
with a = 1 at r ^ Tc. A power-law region develops below the radius 
where Aa = (near rc), with a slope that grows slowly from unity to 
the asymptotic value. 



scenario is based on enhanced feedback due to tidal compression. Consider a 
satellite made of DM and ~ 10% baryons passing through the halo core towards 
a turn-around on the other side. Assume that the baryons have already cooled 
and contracted into the satellite center. The satellite loses its outer DM layers 
in the outer halo such that when it enters the halo core it is baryon-rich. In 
the core, the tides compress the satellite, creating shocks and an efficient burst 
of star formation. Before turn-around on the other side, there is time for the 
resulting supernovae to blow out the satellite gas. If the remaining satellite 
loses half its mass in this blow-out, its density drops by a factor between 8 to 
infinity, depending on whether the gas expulsion is adiabatic or impulsive. Thus, 
the remaining satellite becomes much more susceptible to tidal stripping, which 
disrupts it completely before it re-enters the halo core. 

In our recent work, we address the different problems within the successful 
cosmological framework of CDM by appealing to inevitable feedback effects. In 
Mailer & Dekel (2002) we address the angular-momentum catastrophe, where 
simulations including gas produce disks significantly smaller than the galactic 
disks observed, and with a different internal distribution of angular momentum. 
We first construct a toy model for the angular-momentum buildup by mergers 
based on tidal stripping and dynamical friction, which helps us understand the 
origin of the spin problem as a result of over-cooling in satellites. We then incor- 
porate a simple model of feedback, motivated by Dekel & Silk (1986), and find 
that it can remedy the discrepancies, and in particular explain the low baryon 
fraction and angular-momentum profiles in dwarf disk galaxies. Feedback effects 
may also provide the cure to the missing dwarf problem, where the predicted 
number of dwarf halos in CDM is much larger than the observed number of 
dwarf galaxies (Bullock, Kravtsov & Weinberg 2000). 
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Another approach (e.g., Hogan & Dalcanton 2000) is to appeal to a Warm 
Dark Matter (WDM) scenario, despite the fact that it requires fine-tuning of 
the particle mass to ~ 1 keV. The main feature of WDM is the partial sup- 
pression of small halos, which should help maintaining a core, but simulations of 
halos in WDM seem to still show inner cusps (Bullock, Kravtsov & Colin 2001). 
While the explicit merger picture modeled above may be invalid in this case, 
the gravitational processes involved in the halo buildup still mimic a similar 
behavior. We note that the tidal compression in the core may amplify density 
perturbations and make them behave like merging satellites. 

The suppression of small halos in WDM may remedy the spin catastrophe 
by avoiding over-cooling (Sommer-Larsen & Dolgov 2001), but it is harder to see 
how it would explain the angular-momentum profile in galaxies. Furthermore, 
while the number of dwarfs is already suppressed in WDM, the addition of 
minimum feedback effects is likely to cause an overkill, where the number of 
dwarfs is predicted to be much lower than observed (J. Bullock, private comm.). 

The success of our toy model in matching several independent observations 
indicates that it indeed captures the relevant elements of the complex processes 
involved, and in particular that feedback effects may indeed provide the cure to 
all three problems of galaxy formation in CDM. 
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